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The transition from unicellular to multicellular life was one of a few major events in the history 

of life that created new opportunities for more complex biological systems to evolve. Predation 

is hypothesized as one selective pressure that may have driven the evolution of multicellularity. 

Here we show that de novo origins of simple multicellularity can evolve in response to predation. 

We subjected outcrossed populations of the unicellular green alga Chlamydomonas reinhardtii to 
selection by the filter-feeding predator Paramecium tetraurelia. Two of five experimental populations 
evolved multicellular structures not observed in unselected control populations within ~750 asexual 
generations. Considerable variation exists in the evolved multicellular life cycles, with both cell number 
and propagule size varying among isolates. Survival assays show that evolved multicellular traits 
provide effective protection against predation. These results support the hypothesis that selection 
imposed by predators may have played a role in some origins of multicellularity. 


Nearly all macroscopic life is multicellular. As Leo Buss emphasized in The Evolution of Individuality, the very 
existence of integrated multicellular organisms is an outcome of evolutionary processes, not a starting condition’. 
It seems, in fact, to be a common outcome: multicellular organisms have evolved from unicellular ancestors 
dozens of times*~*. Animals, land plants, fungi, red algae, brown algae, several groups of green algae, cellular and 
acrasid slime molds, and colonial ciliates, among others, each descend from a different unicellular ancestor*”. 

Two retrospective approaches, comparative methods and the fossil record, have proven valuable in recon- 
structing how these transitions may have occurred. Although both approaches have been critical to our under- 
standing of early multicellular evolution, each has its limitations. For most multicellular groups, little or no fossil 
evidence exists that is relevant to the first steps in the transition from unicellular to multicellular life. Comparative 
methods suffer from a lack of intermediate forms between the multicellular organisms we are interested in and 
their extant unicellular relatives. Furthermore, direct knowledge of unicellular ancestors is not available. Extant 
unicellular relatives often serve as stand-ins, but this is a poor approximation, as they have been evolving inde- 
pendently as single-celled organisms since they diverged from their multicellular relatives. 

A third, prospective, approach designed to circumvent these limitations has emerged over the past decade. 
The experimental evolution of multicellularity in otherwise unicellular microbes enables real-time observations 
of morphological, developmental, and genetic changes that attend the transition to multicellular life. Boraas and 
colleagues exposed cultures of the green alga Chlorella vulgaris to predation by the flagellate Ochromonas valles- 
cia, resulting in the evolution of small, heritably stable algal colonies®. Becks and colleagues showed that exposure 
to the predatory rotifer Brachionus calyciflorus selected for heritable changes in the rate of formation of multicel- 
lular palmelloids in the green alga Chlamydomonas reinhardtii’. Ratcliff and colleagues have shown that selection 
for an increased rate of settling out of liquid suspension consistently results in the evolution of multicellular 
‘snowflake’ colonies in the yeast Saccharomyces cerevisiae® and also results in the evolution of simple multicellular 
structures in C. reinhardtii’. 

Predation has long been hypothesized as an ecologically plausible selective pressure that could explain at least 
some origins of multicellularity, as most predators can only consume prey within a narrow range of sizes*'°""’. 
Filter-feeding predators are common in aquatic ecosystems, and algae that are larger than a threshold size are 
largely immune to them”. David Kirk proposed that this ‘predation threshold’ could explain the evolution of 
multicellularity in the volvocine algae, which are close relatives of C. reinhardtii’. 
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Figure 1. Scanning electron micrographs of representative multicellular colonies from evolved populations. 
(A) Shows an amorphous cluster from population B2. Cell number varies greatly between clusters in this clone 
and between clones in this population. (B) Shows an eight-celled cluster from population B5. Octads were 
frequently observed in both populations. 
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In this study, we present experiments in which we used the ciliate predator Paramecium tetraurelia to select 
for the de novo evolution of multicellularity in outcrossed populations of C. reinhardtii. We describe the heritable 
multicellular life cycles that evolved and compare them to the ancestral, unicellular life cycle. In addition, we 
show that the evolved multicellular life cycles are stable over thousands of asexual generations in the absence of 
predators. Comparative assays show that the evolved multicellular phenotypes provide a fitness advantage over 
unicellular algae in the presence of predators. Because C. reinhardtii has no multicellular ancestors, these experi- 
ments represent a completely novel origin of obligate multicellularity’*”. 


Results 

Experimental Evolution. After 50 weekly transfers (~750 generations), simple multicellular structures 
evolved in two of five predator-selected populations (B2 and B5). Such multicellular structures were not observed 
in any of the control populations. Eight strains were isolated from each of three populations (B2, B5, K1). We 
focused our analyses on five focal strains from B2 (B2-01, B2-03, B2-04, B2-10, B2-11) and two strains from B5 
(B5-05, B5-06). Of the isolates from the control population that evolved in the absence of predators (K1), we 
analyzed two strains (K1-01, K1-06). Phenotypes of other isolates from populations B2, B5 and K1 did not differ 
qualitatively from the focal strains and were not investigated further. The strains have maintained their evolved 
characteristics of simple multicellularity in the absence of predators for four years as unfrozen, in-use laboratory 
strains. Therefore, we are confident that the phenotypic traits that we report below are stably heritable. 

Some strains, notably those from population B2, appeared to form amorphous clusters of variable cell number 
(Fig. 1A). Other strains, notably those from population B5, commonly formed stereotypic eight-celled clusters, 
with an apparent unicellular and tetrad life stage (Fig. 1B). Other phenotypic differences could be easily discerned 
by light microscopy. For example, in Fig. 1, an external membrane is visible around both evolved multicellu- 
lar colonies, indicating that they formed clonally via repeated cell division within the cluster, rather than via 
aggregation. 

We found a variety of life cycles within and among populations. Using time-lapse videos of evolved strains 
growing in liquid culture in 96-well tissue culture plates, we qualitatively classified strains into four life cycle 
categories to compare their similarities and differences (Fig. 2A—D). Sorting observations of a complex trait such 
as a life cycle into categories unavoidably requires some subjectivity, but our category assignments are in good 
agreement with quantitative data presented in the next two subsections. We discuss each life cycle category below, 
and then present quantitative measurements on parent and propagule cluster sizes. 

In population K1, which evolved without predators, ancestral life cycle characteristics of the unicellular, 
wild-type C. reinhardtii were retained (Fig. 2A, Supplemental Videos 1 and 2). Specifically, as cells reproduce 
asexually, they lose motility and undergo 2-5 rounds of mitosis before releasing motile, single-celled propagules. 
It should be noted that even in wild-type C reinhardtii, the dividing parent cluster is a transient multicellular stage; 
however, it does not persist after propagules are released. Interestingly, in the two populations that evolved multi- 
cellularity in response to predation (B2, B5), strains B5-06 and B2-10 retained a life cycle typical of the ancestral, 
wild-type C. reinhardtii (Supplemental Videos 3 and 4, respectively). 

Life cycles of the remaining strains isolated from populations B2 and B5 were distinct from wild-type 
Chlamydomonas, as clusters of various sizes persist through multiple rounds of reproduction. Ordinarily, strain 
B2-01 releases motile, single-celled propagules during reproduction, similar to wild-type (Supplemental Video 5). 
However, in some clusters, cells undergoing division separate, but remain proximately located because they are 
embedded in an extracellular matrix (ECM) of the parent cluster (Fig. 2B). As these cells continue to grow and 
divide, some remain embedded in the ECM, which creates growing aggregations of cells. Strain B5-05 also pro- 
duces motile, single-celled propagules that are often embedded in the maternal ECM (Supplemental Video 6). In 
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Figure 2. Depiction of C. reinhardtii life cycles following evolution with (B2, B5) or without (K1) predators for 
50 weeks. Categories (A-D) show a variety of life cycle characteristics, from unicellular to various multicellular 
forms. Briefly, A shows the ancestral, wild-type life cycle; in B this is modified with cells embedded in an 
extracellular matrix; C is similar to B but forms much larger multicellular structures; while D shows a fully 
multicellular life cycle in which multicellular clusters release multicellular propagules. Evolved strains were 
qualitatively categorized based on growth during 72-hour time-lapse videos. Strains within each life cycle 
category are listed below illustrations. Representative microscopic images of each life cycle category are at the 
bottom (Depicted strain in boldface). 


addition to retaining propagules embedded in the ECM, growing clusters of B5-05 ensnare free-swimming cells, 
creating aggregations that grow much larger than those of B2-01 (Fig. 2C). 

Conversely, three of the strains isolated from population B2 exist in cell clusters comprised only of direct 
descendants, as opposed to chimeric aggregations with free-swimming cells. Clusters from strains B2-11, B2-03, 
and B2-04 grow in tightly associated groups of direct descendants embedded in the maternal cell wall (Fig. 2D; 
Supplemental Videos 7, 8 and 9, respectively). Development in these isolates is therefore strictly clonal, with 
important implications for evolvability. Since the cells within a multicellular structure are likely to be genetically 
identical, other than differences resulting from new mutations, genetic variation in a population would be parti- 
tioned primarily among colonies. The clonal development observed in these isolates therefore suggests that the 
observed multicellular clusters would be well-suited to serving as units of selection. 

In key respects, the isolates from population B2 appear to have recapitulated early steps hypothesized as 
leading to differentiated multicellularity in volvocine algae’®. In fact, the evolved multicellular algae are sim- 
ilar in their gross morphology to small colonial volvocine algae such as Pandorina. Furthermore, a degree of 
genetic control of cell number similar to that seen in undifferentiated colonial volvocine algae manifests in our 
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Figure 3. Median cluster sizes of evolved strains at the time-point where strains reached maximum size. To 
determine the time where strains were largest, we calculated the mean from five replicate medians (open dots), 
for each strain and time-point. Because data are not normally distributed, medians were chosen to approximate 
central tendency. Cells per cluster were measured by sampling strains over six days of growth, staining nuclei 
with DAPI, and imaging using fluorescent microscopy. From left to right, time-points of maximum size for each 
strain were: 12, 12, 72, 120, 108, 84, 96, 96, and 72 hours. From left to right, median sizes of replicate populations 
averaged 1, 1, 1.4, 1.8, 1.8, 4, 6.2, 5.6 and 5.6 cells per cluster. Letters above columns indicate significant 
differences among medians (Kruskal-Wallis, X; = 38.49, p <0.001, Dunn’s Test for multiple comparisons, 

a= 0.05). Sizes at the initial time-point (0 hrs) were omitted from analysis because they represent starting 
conditions. Shading is only for ease of visualization. 


eoceeeree eee eee eee eee eee eee eee eee eee ee ee eee eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeseeeseeeeeeeeeseeeeeeseeeeeeeeeeeeeeeeeeeeeeeeeeseeeeeeeeeeeeeoee 


evolved multicellular strains. Propagules of these strains are typically multicellular and, critically, no motile prop- 
agules were observed in these strains. Propagules from the evolved multicellular strains were nearly all 4-, 8-, or 
16-celled, a range similar to that of small colonial volvocine algae and smaller than that of Pleodorina starrii, in 
which propagules of a single genotype can span a 16-fold range of cell numbers’’. 


Maximum Cells Per Cluster. In order to determine the size attained by evolved strains, we sampled, stained, 
and imaged growing cultures in tubes every 12 hours over a 6-day period. We found that cluster sizes of strains 
that evolved without predators (K1-01, K1-06) and strains that retained ancestral life cycle traits, despite evolving 
with predators (B5-06, B2-10), remain small throughout the six days of growth (Fig. 3). Even at their largest, 
median cluster sizes of replicate populations of strains K1-01 and K1-06 were all one cell per cluster, showing that 
these strains exist as unicellular individuals with a transient a multicellular stage occurring during reproduction. 
Similarly, we find that cells of strain B2-01 exist primarily as single cells, with median sizes averaging 1.8 cells. 
Under time-lapse observation, we found that dividing cells of B2-01 would occasionally separate from direct 
contact with each other while still being held together by maternal ECM/cell wall. This was seen as white spaces 
appearing between cells, which then remained equidistant in a rigid structure. During the procedures to stain 
and image cells, it is possible that these partly separated cells were dispersed, leading to a slightly lower estimate 
of median cluster sizes. The median cluster sizes of replicate populations of strain B5-05 averaged 4.0 cells, while 
those that retained the ancestral phenotype (K1-01, K1-06, B5-06, and B2-10) averaged 1.0, 1.0, 1.4, and 1.8 cells, 
respectively, and B2-01 averaged 1.8 cells. Unlike strains B2-01 and B5-05, in which ECM holds cells in clusters, 
cells of the largest strains are visibly encapsulated within the mother cell wall. Medians of replicate populations of 
these strains (B2-11, B2-03, and B2-04) averaged 6.2, 5.6, and 5.6, respectively. 

Because cells were mixed several times during the staining and imaging process, it is possible that some clus- 
ters of cells were disrupted. Thus, cluster sizes deduced by fluorescence microscopy are likely to be underesti- 
mates, possibly affecting cluster size estimates for certain strains more than others (i.e., those held together by 
extracellular matrix rather than by maternal cell wall). 


Number of Cells per Propagule. In order to determine propagule sizes of evolved strains, we manually 
analyzed 72-hour time-lapse videos, recording numbers of cells in propagules released from parent clusters. 

For most strains the majority of propagules were single-celled, as observations were skewed toward smaller 
propagules, which require less biomass to produce (Fig. 4). To accurately depict the context of individual cells in 
propagules, we also show the data as boxplots, weighted by biomass (Fig. 4, black boxplots). Strains that main- 
tained ancestral characteristics (K1-01, K1-06, B5-06, B2-10) almost exclusively released single-celled propagules 
(zero out of N> 150 propagules were multicellular; Fig. 4). Although cells from strains B2-01 and B5-05 often 
exist in clusters, the vast majority of propagules released by these reproducing clusters were single-celled as well. 
In strains with larger average maximum cluster sizes (B2-11, B2-03, B2-04), cells are more frequently released in 
multicellular propagules of up to 64 cells in size. Thus, among the strains that evolved under predation, a variety 
of size-related traits emerged in both cluster size and propagule size. 
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Figure 4. Sizes of propagules released by evolved strains during a 72-hour time-lapse. Propagule sizes were 
manually measured from time-lapse videos. Green dots show actual observations of propagule size. Because the 
frequency of observations is skewed toward smaller propagule sizes, we show the data as weighted boxplots as 
well, where a 32-celled propagule has 32 times the weight of a single-celled propagule. Letters above columns 
indicate significant differences among medians (Kruskal-Wallis, X; = 2779, p < 0.001, Dunn's Test for multiple 
comparisons, ~a = 0.05). Sample sizes of propagule observations for each strain are indicated along the bottom of 
the plot. Strains are ordered by weighted-mean propagule size. 
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Figure 5. Average absorbance value differences (A,,,,) between populations with and without predation 

for evolved multicellular and unicellular C. reinhardtii strains. Multicellular strains averaged a much lower 
absorbance value difference over the duration of the experiment than unicellular strains. Error bars are standard 
deviations of twelve strains per treatment. Markers are offset 0.5 h so that error bars can be seen. 
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Defense Against Predators. Multicellularity appears to provide an effective defense against predation by 
the rotifer Brachionus calyciflorus (our choice of predator for this assay is explained in the Methods). Comparisons 
of predation rates were made by measuring the absorbance of each of the four replicates without predators and 
subtracting the average of the four replicates with predators. This value, labeled A,,,, shows markedly different 
trajectories for unicellular and multicellular treatments with predators (Fig. 5). To simplify the analysis, infer- 
ences were drawn from the final time-point measurements (hour 117). A test for an effect of phenotype (unicellu- 
lar or multicellular) on A,,, values at hour 117 shows that phenotype has a strong effect on algal density (ANOVA, 
F, 95 = 7.37, P=0.013). The unicellular strains, on average, experienced ~2.5-fold greater rates of predation com- 
pared to the multicellular strains. That is, the mean A,,, at hour 117 (day 5) were 0.367 (s.d. 0.125) and 0.140 (s.d. 
0.260) for 12 unicellular and 12 multicellular strains, respectively. Over the time course, the slopes of A,,, dif- 
fered significantly between unicellular (0.0030) and multicellular (0.0002) strains using either linear or quadratic 
regression (P < 0.0001 for the interaction between phenotype and time in both analyses). 





SCIENTIFIC REPORTS | (2019) 9:2328 | https://doi.org/10.1038/s41598-019-39558-8 5 


www.nature.com/scientificreports/ 





Discussion 

Our results show that the transition to a simple multicellular life cycle can happen rapidly in response to an eco- 
logically relevant selective pressure. By increasing in size beyond the “predation threshold”? of a filter-feeding 
predator, multicellular C. reinhardtii that evolved from an ancestrally unicellular lineage are protected from pre- 
dation for at least part of their life cycle. Under selection for increased size, formation of multicellular structures 
may be an easier route than increasing cell size because of trade-offs imposed by scaling relationships (chiefly 
the reduction in surface-area-to-volume ratio)'®, because more mutational paths are available, and/or because 
available mutations have fewer or less severe pleiotropic effects. 

Observations in other species suggest that the relative ease of transitioning from a unicellular to a multicellu- 
lar life cycle is at least somewhat general. Similar transitions reportedly occurred within 100 generations in the 
green alga Chlorella vulgaris® and within 300 generations in the yeast Saccharomyces cerevisiae®. Although there 
have surely been many microbial evolution experiments in which multicellularity did not evolve, we can at least 
be confident that this phenomenon is not unique to Chlamydomonas. Grosberg and Strathmann‘ may have been 
right to call the evolution of multicellularity a “minor major transition’. 

However, the available evidence also suggests a substantial stochastic component to the evolution of mul- 
ticellularity. In previous experiments using settling selection as opposed to predation selection, multicellular 
structures evolved in one of ten selected populations in C. reinhardtii’ and in “about 70%” of “many” selected 
populations in C. vulgaris®. In the experiment reported here, a variety of multicellular forms evolved in two of 
five selected populations. Only in S. cerevisiae has the evolution of such forms proven consistent across replicate 
populations®. 

The multicellular forms that evolved in our experiments bear some resemblance to palmelloids, transient mul- 
ticellular structures that appear under some conditions in C. reinhardtii and related algae. The propensity to form 
palmelloids varies among species and strains’®”’. In wild-type C. reinhardtii, it is a phenotypically plastic response 
to the presence of certain organic acids*', Ca** deficiency, chelating agents, high phosphate concentrations”, and 
certain zooplankton predators***. Predators have also been observed to induce the formation of multicellular 
structures in the green algae Chlorella*”° and Scenedesmus*’”*, as well as in bacteria”””®. 

The term “palmelloid” encompasses two distinct phenomena: a failure of daughter cells to escape the mother 
cell wall, resulting in ‘colonies’ of four to sixteen cells, and the production of a gelatinous material that is secreted 
extracellularly, resulting in adhesion among cells and/or colonies”°. Cell walls are easily distinguished from ECM 
under magnification, since cell walls block transmitted light much more strongly (i.e. appear much darker) than 
ECM. We observe both phenomena in evolved multicellular isolates. 

The ability of wild-type C. reinhardtii to form palmelloids suggests that the founding population in our exper- 
iment already possessed a toolkit for producing multicellular structures. However, while the palmelloid condition 
is expressed facultatively in wild-type C. reinhardtii, the strains that evolved in our experiment are obligately mul- 
ticellular. Like palmelloids, our evolved multicellular isolates lack motility, suggesting that the ability to express 
both unicellular and multicellular phenotypes may be optimal when predation pressure varies over time”’. 

Our life cycle observations, carried out in as near as possible identical conditions, constitute a common-garden 
experiment, demonstrating that the evolved phenotypes have a genetic basis. It may be, though, that this basis 
involves the co-option of a previously existing plastic response. If so, the shift from a primarily unicellular (but 
facultatively multicellular) to an obligately multicellular life cycle may have required only a change from faculta- 
tive to obligate expression of the genes involved in palmelloid formation. 

Our results support the view that predation may have played a role in at least some origins of multicellularity. 
Consistent with prior results®’, predation can drive the evolution of simple multicellular structures, and these 
structures provide protection from predators®**. The ‘chicken-and-egg’ problem brought up by some authors 
— that multicellular predators are required to drive the evolution of multicellularity”'* — does not exist: unicel- 
lular filter-feeding protists exist, including, for example, Paramecium, and there is no reason to think that such 
predators only arose after animals. Furthermore, some origins of multicellularity, including those of brown and 
volvocine algae, postdate that of animals, and could therefore have been driven by animal predation. 


Methods 

Experimental Evolution. To increase genetic diversity at the onset of the experiment, we founded five 
experimental (B1-B5) and three control (K1-K3) populations from an F1 cross. We obtained plus (CC-1690) 
and minus (CC-2290) mating type strains of C. reinhardtii from the Chlamydomonas Resource Center (www. 
chlamycollection.org) and crossed them using a protocol adapted from Harris’’. Vegetative cells were grown to 
high density in 10 mL TAP medium”? at 22.5°C and pelleted in 15 mL conical tubes by centrifugation for 5 min. 
at 2000 x g. Pelleted cells were resuspended in 10 mL mating medium (TAP minus NH,Cl) and grown in light 
for 6h to induce gametogenesis. Gametes were mixed in a 24-well tissue culture plate with 0.5 mL of each mat- 
ing type, incubated in light for 4h, and then dried in the dark. Single mating-type controls of 1 mL were treated 
identically to ensure that no vegetative material survived desiccation. After 20 days in the dark, crosses and single 
mating-type controls were each flooded with 1 mL TAP medium and placed in light at 22.5 °C to germinate. Single 
mating-type controls showed no evidence of vegetative growth after two weeks; thus, desiccation was effective in 
killing vegetative material and surviving algae from the crosses must be outcrossed. 

Experimental populations were started by mixing 0.1 mL of C. reinhardtii F1 cultures, grown to ~2 x 10° cells/ 
mL, with 0.1 mL ofa culture of P. tetraurelia, grown to ~2 x 10° cells/mL, in 1.3 mL of COMBO medium”, mak- 
ing a final volume of 1.5 mL. Control populations were founded by inoculating 0.1 mL of C. reinhardtii F1 cultures 
in 1.4mL of COMBO medium. All populations were propagated in 24-well plates and transferred weekly. 

During transfers, populations were homogenized with a multichannel pipette by drawing and dispensing 
1 mL three times. After mixing, 0.1 mL of culture was transferred to 1.4 mL fresh COMBO in a new 24-well plate. 
After 350 days, populations were plated on TAP medium +1.5% agar. For each population plate, eight individual 
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colonies were randomly picked and re-plated on agar. Individual colonies were picked and re-plated at least three 
times each to ensure that each was monoclonal and free of contaminants. 

We initially characterized evolved phenotypes using light microscopy. A homogenized dilution of each isolate 
grown in TAP medium was fixed in iodine and stored at 4°C. For each isolate, 48 uL samples were photographed 
in replicate on a hemocytometer. Images for each of the clones were processed in ImageJ to determine particle 
counts and relative particle sizes, and these data were later processed in Python to determine the relative pro- 
portions of different sized particles. The size threshold used to distinguish between unicellular and multicellular 
particles was user-determined and was rarely violated by any samples. In all instances where the size threshold 
was violated, small multicellular particles were categorized as unicellular; thus, we remain conservative in our 
characterization of particle sizes, under- rather than overestimating the number of cells in particles. 


Maximum Cells Per Cluster. In order to determine the size attained by different evolved strains, measured 
as the maximum number of cells per cluster, we stained and imaged samples via fluorescent microscopy over a 
6-day examination period. Strains were cultured in 10 mL of TAP medium in a climate-controlled chamber on a 
14:10-hour light:dark cycle (23.5°C, 41% humidity, 326 lx). In order to standardize starting cell concentrations, 
preliminary cultures were first inoculated from scrapes of colonies on TAP plates (2% agar) and allowed to grow, 
without disruption, for 6 days until reaching stationary phase. Cultures were then transferred 1:100 into fresh 
TAP medium and examined over the following 6-day examination period. 

Throughout the examination period, cultures were sampled, stained, and imaged every 12 hours in order to 
capture the maximum size attained by various strains before reaching stationary phase. Strains B5-05 and B5-06 
were not sampled at the final time-point. Samples were extracted from well-mixed culture tubes and transferred 
into 1.5mL microcentrifuge tubes (200 1L from the first 6 time-points to compensate for low cell density, 100 1L 
from all time-points thereafter). The samples were centrifuged at 250 x g for 5 minutes and the supernatant was 
replaced by a 50% v/v ethanol-water solution to fix cells and make cell membranes more permeable to the DAPI 
(4’,6-diamidino-2-phenylindole) fluorescent nucleotide stain. After a 10-minute fixation period, samples were 
centrifuged again (250 x g for 5 minutes) and the supernatant was replaced with distilled water containing 1 g/ 
mL DAPI. Cells with DAPI were kept in darkness for one hour and then moved to a refrigerator without light to 
prevent photo bleaching until imaging on the microscope. 

Prior to imaging, samples were removed from the refrigerator and inverted to mix cultures. Seven L of each 
sample were mounted on glass microscope slides and large 6 x 6 images were stitched from individual images 
captured at 40x magnification, creating high-resolution (441 megapixel) wide-field images, providing a robust 
measurement of the distribution of clusters sizes within the population. 


Image and Data Analysis. Large stitched images of individual samples were imported and analyzed in 
Fiji’. We wrote and implemented a script in ImageJ macro (Supplemental Material) to automatically calculate the 
number of cells in each cluster. The script accomplished this by demarcating cluster boundaries and recording the 
number of pixel maxima within each boundary, corresponding to the number of DAPI-stained nuclei. Following 
this procedure, all images used in our analyses were manually screened and any regions of interest resulting from 
artifacts (e.g., dust particles, autofluorescence) were removed from the dataset. Data for five technical replicates 
of all tested strains, each represented by 12 time-points sampled over six days, were concatenated into one large 
dataset. 

In order to determine the maximum number of cells per cluster achieved by each strain, hereafter referred to 
as “cluster size’, we first determined the time-point at which clusters were largest. This was calculated separately 
for each strain, as different strains reached their maximum size at different time-points during the six-day culture 
cycle. First, we calculated the median cluster size of each replicate at each time-point throughout the culture cycle. 
Medians were chosen to represent the central tendency instead of the mean because the distributions of cluster 
sizes were right-skewed and occasionally contained large outlier clusters. Then, to determine the time-point at 
which cluster size was largest for each strain, we calculated the average (mean) of the median cluster sizes from 
replicate populations within each time-point. The first time-point, representing cluster sizes at time of inocula- 
tion, was removed from analysis. 


Number of Cells Per Propagule. In order to determine the number of cells per propagule released by clusters 
of each strain, hereafter referred to as “propagule size’, we time-lapse imaged all strains over 72 hours of growth. 
First, cells were inoculated from colonies on TAP agar plates into liquid TAP media and incubated for five days 
to grow to high density. Then, each culture was homogenized by vortexing, transferred 1:100 into five replicate 
tubes of fresh TAP media, and allowed to grow for two days. These exponentially growing cultures were mixed 
by vortexing, diluted 1:10 in TAP medium, and then 100\L from each culture was randomly inoculated into the 
central 60 wells of a 96-well tissue culture plate that already contained 100 :L of TAP medium in each well. 

For time-lapse microscopy, the 96-well plate (containing five technical replicates of each strain) was imaged at 
200x magnification, where the field of view was positioned on cells near the center and at the bottom of wells. The 
time-lapse was run for 72 hours, capturing images of each well every 30 minutes. 

Time-lapse images were manually analyzed using Image]. For each well, all images from the time-lapse were 
viewed individually and sequentially. All cells or clusters of cells in the frame at the beginning of the time-lapse 
were labeled with a number indicating their identity. Throughout the time-lapse, when a cell or cluster of cells was 
seen to separate from the initial parental cluster/group, the number of cells in the propagule was recorded. For 
wells that had few or no initial cells in the frame (typically occurring for motile, unicellular strains), cells or clus- 
ters of cells were labeled as they entered the field of view, until collecting a sufficient sample size (N > 5 parental 
clusters). Due to the long time required for clusters to fully separate physically—even after they have detached 
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from a parent cluster—propagule sizes were recorded four frames (two hours) after the propagule was observed 
to initially split from the parental cluster. Data collected from time-lapse images were analyzed in R*®. 


Defense Against Predators. To test whether multicellularity affords protection from predators, we sub- 
jected evolved unicellular and multicellular isolates to predation by the rotifer Brachionus calyciflorus. Rotifers 
were chosen for this assay instead of Paramecium to reduce the likelihood that some adaptation other than cluster 
formation protected cluster-forming strains from the predator. Rotifers were obtained as cysts in vials and stored 
at 4°C in the dark. To ensure that cultures of B. calyciflorus were free of contaminants, cysts were treated and 
incubated according to the protocol described by Suga et al.°’. Neonate rotifers were grown axenically in 10 mL 
of WC medium”* with 100,.g/mL ampicillin in culture tubes at the same temperature and light conditions as the 
C. reinhardtii cultures. In order to sustain the population of B. calyciflorus before the experiment, 1.0 mL of an 
axenic, unicellular strain of C. reinhardtii, CC-1690, was inoculated into the culture via micropipette. 

Evolved isolates of C. reinhardtii were grown to high density, and culture absorbance was measured using a 
microplate reader (Molecular Devices SpectraMax® M5, 420 nm), whereafter each was diluted to approximately 
3.45 x 10° cells/mL. Four replicates of each strain were randomly assigned to 24-well plates and 1.33 mL of each 
strain pipetted into designated wells. Twelve multicellular strains (B2-03, 04, 06, 07, and 11; B5-01 through 05, 07, 
and 08) and twelve unicellular strains (B2-01, 10, and 12; B5-06; K1-01, 02, and 04 through 09) were used for this 
assay. Calibration curves relating absorbance to cell density were established for each strain by comparing absorb- 
ance readings to hemocytometer counts. Because our experimental design compared differences of absorbance 
readings with and without predators for each strain, it is robust to among strain differences in absorbance per cell. 
That is, it does not require that equal absorbance readings indicate equal cell densities across strains, only that 
absorbance readings scale linearly with cell density within each strain. 

Rotifer cultures were poured through a 251m mesh filter to remove the algae they had been feeding on, then 
pipetted into a sterile test tube awaiting transfer to the experimental wells (C. reinhardtii cells are ~10 1m in diam- 
eter and easily pass through the filter). Depending on the treatment condition, 0.67 mL of either WC medium 
or WC medium with predators was added to each well. Initial absorbance values were recorded and then plates 
were stored under lights in an incubator on a 14/10-hour light/dark cycle at 22.5°C for 5 days. Throughout the 
14-hour light period, absorbance readings were taken every 3 hours, providing growth rates with high resolution 
for all strains. 


Data Availability 
The datasets generated during and/or analysed during the current study are available in the Dryd repository, 
https://doi.org/10.5061/dryad.53k345s. 
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